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Abstract. Interactions between dark matter and dark energy which result in a power-law 
behavior (with respect to the cosmic scale factor) of the ratio between the energy densities 
of the dark components (thus generalizing the ACDM model) have been considered as an 
attempt to alleviate the cosmic coincidence problem phenomenologically. We generalize this 
approach by allowing for a variable equation of state for the dark energy within the CPL- 
parametrization. Based on analytic solutions for the Hubble rate and using the Constitution 
and Union2 SNIa sets, we present a statistical analysis and classify different interacting and 
non-interacting models according to the Akaike (AIC) and the Bayesian (BIC) information 
criteria. We do not find noticeable evidence for an alleviation of the coincidence problem 
with the mentioned type of interaction. 

Keywords: Dark Matter, Dark Energy, Interacting models, CPL parametrization. 



Contents 



1 

_L 


Introduction 


i 

J. 


2 


A suitable interaction 


2 




2.1 


General relations 


2 




2.2 


The non-interacting case 


4 




2.3 


The case £ = 1 


4 




2.4 


The case £ = 3 


4 


3 


Statistical analysis and observational constraints 


6 




3.1 


Observational tests 









3.1.1 SNIa 


6 






3.1.2 BAO 


7 






3.1.3 CMB 


8 




3.2 


Analysis of the models 


8 






3.2.1 The ACDM model 


8 






3.2.2 The non-interacting case 


9 






3.2.3 The case £ = 1 


11 






3.2.4 The case £ = 3 


11 






3.2.5 Comparing the models 


13 


4 


Summary 


16 



1 Introduction 

Despite of the many efforts of the past decade, the physical nature of the cosmological dark 
sector remains mysterious. The preferred model for the current cosmological dynamics is 
still the ACDM scenario, although there is an ongoing discussion (see, e.g., [1]) on shortcom- 
ings of this model which, however, became a reference model also for competing approaches. 
Much work has been done to study deviations from the ACDM model. Usually, this implies 
introducing new parameters which are then constrained by the ever increasing number of 
available observational data. In this way one tries to clarify, e.g., whether the effective dark- 
energy equation-of-state (EoS) is really described by a constant parameter and whether this 
parameter is really equal to -1 as required by a true cosmological constant. Most investiga- 
tions impose limits on potential deviations from ACDM, but generally the latter turns out 
to be consistent with the data and even remains the preferred option. For recent results see, 
e.g., [2] and [3]. 

The ACDM model is characterized by a pressureless dark matter component with an 
energy density p m which decays with the third power a 3 of the cosmic scale factor a and 
a constant energy density p\, associated with the cosmological constant. Consequently, for 
the ratio of the energy densities one has ^ oc a~ 3 . Given the huge range of values for a 
during the cosmic evolution, a currently measured ratio ^ = 0(1) seems to single out the 
present epoch as a very special period. This is known as the cosmic coincidence problem. 
There have been several attempts which tried to uncover a dynamical mechanism to make 
this coincidence in a sense natural. One of the possibilities that have been studied over the 
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last years is to admit an interaction between dark matter and dark energy and to check its 
influence on the ratio of the energy densities of both components. In this paper we consider 
an approach that relies on a phenomenological relation for the ratio of the energy densities 
of dark matter, p m , and a dynamical dark energy component with density p x , 

Pm —£ /-i -i \ 

— = ra ? , (1.1) 

Px 

where the scale factor a of the Robertson- Walker metric was normalized to a present value 
a(to) = cio = 1, and r is the energy-density ratio at the present time. The power £ is a 
constant, non-negative parameter. A behavior of the type (1.1) was suggested by Dalai et 
al. [4] in order to address the coincidence problem on a phenomenological basis. For an 
equation of state p x = —p x of the dark-energy component, a value £ = 3 amounts to the 
ACDM model. A value £ = represents a stationary ratio y 1 = const. According to [4], the 
deviation of the parameter £ from £ = quantifies the severity of the coincidence problem. 
Any solution with a scaling parameter £ < 3 will make the coincidence problem less severe. 
Investigations along this line include [5-7] . Here we generalize previous studies by admitting a 
time varying equation of state for the dark-energy component. Using the widely applied CPL 
parametrization [8] for the equation of state parameter, we obtain analytic solutions for the 
special cases £ = 1 and £ = 3. In a subsequent statistical analysis our results are confronted 
with the SNIa data from the Constitution [9] and Union2 [10] samples and complementary 
information from the baryon acoustic oscillations (BAO) scale [12] and the position of the 
first acoustic peak in the CMB TT spectrum [13]. Similarities and differences to alternative 
studies in the literature are pointed out. The results are assessed according to the Akaike 
(AIC) and Bayesian {BIC) information criteria. We find that none of the models considered 
here can really compete with the ACDM model. However, this conclusion is partially based 
on a prior which is adapted to the ACDM model so that the latter seems naturally preferred. 
Therefore, the mentioned result may not yet be the final answer. But at least at the moment, 
there is no evidence that the coincidence problem can be alleviated by an approach based on 
(1.1) with a corresponding interaction. 

The paper is organized as follows. In section 2 we present our interacting dark-energy 
model and obtain analytic solutions for the Hubble rates for two special cases. A correspond- 
ing statistical analysis is the subject of section 3. A summary and a discussion of our results 
are given in section 4. 

2 A suitable interaction 
2.1 General relations 

Let us assume the cosmic medium to be dynamically dominated by cold dark matter (subindex 
m) and dark energy (subindex x). Additionally, we take into account a baryon component 
(subindex b). Then 

P = Pm + Px + Pb and p = p m +Px+Pb (2.1) 

are the total energy density and the total pressure, respectively. The components are assumed 
to possess the equations of state 

P™« Pm , Pb<^ Pb and p x = wp x . (2.2) 



- 2 - 



An interaction between both dark components may be described by the set 



Pm + 3Hp m = Q 



(2.3) 



and 

p x + 3H (l + u)p x = -Q . (2.4) 
For the separately conserved baryon component we have 

p b + 3Hp b = => p b oc a" 3 . (2.5) 

A specific expression for the interaction term can be derived by considering the time evolution 
of the ratio — , 



Pm 
Px 



Pm 


Pm 


Px 


= 3H^ 


UJ + 


Px 


Pm 


Px_ 


Px 





Pm + Px Q 
PmPx 3H 



(2.6) 



Inserting here the ansatz (1.1) and solving for the interaction term yields 



Q 



00 + 



3Hp, n 



ra 



+ 1 



(2.7) 



Eq. (2.7) demonstrates that by choosing a suitable interaction between both components, we 
may produce any desired scaling behavior of the energy densities [6]. The uncoupled case, 
corresponding to Q = 0, is given by | + uj = 0. The ACDM model corresponds to w = — 1 
and £ = 3. Generally, interacting models are parameterized by deviations from £ = — 3w. 
Any solution which deviates from £ = —3w represents a testable, non-standard cosmological 
model. For £ > the expression (2.7) becomes very small for a <C 1. Consequently, the 
interaction is not relevant at high redshifts. This guarantees the existence of an early matter- 
dominated epoch. An energy transfer from dark energy to dark matter, i.e. Q > 0, requires 
§ + uj < 0. 

All relations so far are valid for EoS parameters that are not necessarily constant. In 
order to take into account a variable EoS parameter, we resort to the CPL parametrization 
[8] 

lj = ujq + U)i(l — a) , (2-8) 

where ujq and uj\ are constant parameters. Then the condition for Q > depends on a. We 
have ujq + uj\ + | < for a <S 1, ojq + | < for a = 1 and ojq — uj\a + | < for a 3> 1. But 
for the latter limit this kind of parametrization does not seem to be useful. 
The dimensionless interaction quantity (2.7) can be written 



Q 



3Hp, n 



£ ( a^ +1 

, (oj +UJl + -) + Ul 
r + a?/ 3 \r + a« 



(2- 



In the following we consider separately the analytically solvable cases £ = 1 and £ = 3. 
The latter is expected to test primarily deviations from the ACDM model due to a time- 
varying equation of state for the dark energy, the former should provide information about 
the feasibility to alleviate the coincidence problem with the help of a suitable interaction. 

To situate our analysis properly, we start with the non-interacting case in the following 
subsection. 
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2.2 The non-interacting case 

If the dark components are uncoupled, equivalent to Q = in the balances (2.3) and (2.4), 
the resulting square of the Hubble rate, after implementing the CPL parametrization, is [14] 



H(a) 
Hn 



H o/i , . . , , x exp [3wi (a - 1)J , 



where 



3 #o 



a 3(l+wo+^i) 



3flg 



and J) T , 



.t-0 



8vrGp : 



.r0 



3flg 



(2.10) 



(2.11) 



with Q 60 + S^mo + ^xo = 1- 
2.3 The case £ = 1 

Under this condition, integration of the balance equations (2.3), (2.4) and (2.5) provides us 
with 



Pm — PinOd 



r + a 
r + 1 



-3(i+u)o+wi+a;ir) 



exp [3^i (a — 1)] 



Pm 

Px = — a 
r 



and 



Pb = Pboa 



-3 



respectively. Taking into account 



n 



inO 



n 



m.0 



£I x q 1 — f2/,n — f2 



b0 ~~ "mO 



(2.12) 
(2.13) 

(2.14) 



we find 



H(a) 
Hn 



n 



bo 



+ 



(i-n b0 + n m0 zf y a ^+y) 



— exp [3o;i(a-l)] , (2.15) 



for the square of the Hubble rate, where 

y = too + wi 



1-0 



60 



1 — OfeO — ^mO , 

Disregarding the baryon component and setting oj\ = we recover the results of [6]. 
2.4 The case £ = 3 

For this case the analytic solution for the matter density is 



(2.16) 



Pm — PmO a 



r + a" 
r + 1 



3"l -(wo+wi+1) 



^1/3 _|_ -y /j.2/3 — ^1/3q, _j_ ^2 
r l/3 _|_ a y r 2/3 _ r l/3 _|_ i 



,[3ui(o-l)+A] 



where 



A = 3wi 



,1/3 



arctan 



- 1 / 3 - 2a 



arctan 



r i/3 _ 2 
V3rV3 



(2.17) 
(2.18) 
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With the abbreviation u = ojq + oo\ + 1, the corresponding Hubble rate is given by 



H[a) 



l 2 



^feo Q r + a 3 



r + 1 

r + a 3 




V3 + 1 



1/3 



,[3wi(o-l)+A] 



(2.19) 



In the following we shall investigate the cosmological dynamics based on the Hubble rates 
(2.10), (2.15) and (2.19). 

To illustrate the situation with respect to the coincidence problem, we compare the 
behavior of VL m = 8irGp m / '(3H 2 ) and VL X = 87rGp x /{3H 2 ) for the ACDM model, for the 
interacting model studied in Ref. [7] (with a constant EoS parameter u) = —1.01 and £ = 3.16) 
and for our interacting models with £ = 1 and £ = 3 in Fig.l. To plot and to compare the 
corresponding curves qualitatively, we rely on the fiducial best-fit values Q m o = 0.289 and 
uo = —1.01 of Ref. [7], for all the models shown here. For the case £ = 1, both curves are 
considerably closer to each other for a wide redshift range than for any of the other models, 
which we interpret as an alleviation of the coincidence problem. Fig. 1 shows also that for 
£ = 1 the equality of both components occurs earlier in time than for the other models. In 
previous investigations some room was left in the parameter space for this model [5]. We 
reconsider this issue here again on the basis of the recent data from SNIa, from BAO and 
from the CMB shift parameter. 




The 



Figure 1. Redshift dependence of the fractional contributions fl m = 8 3^2'" and Q x — 8 ^§ x 
baryon contribution is neglected. In this illustration the same values fi m o = 0.289 and to = —1.01 
were used for all the models. While the curves for the interacting model of [7] (solid lines) and for 
our £ = 3 model (which is indistinguishable from the ACDM model) are similar to each other, the 
difference between f2 m and £l x is much smaller for £ = 1 than for any of the other models. 
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3 Statistical analysis and observational constraints 



As already mentioned, our aim is to compare models with the Hubble rates (2.10), (2.15) and 
(2.19). To test the different models against cosmological observations, we consider the SNIa 
data sets Constitution [9] and Union2 [10], as well as baryon acoustic oscillations (BAO) [12] 
and the CMB shift parameter [13]. For the SNIa data we adapt the procedure put forward 
in [15]. 

3.1 Observational tests 
3.1.1 SNIa 

Usually, the statistical x 2 analysis procedure is based on the probability distribution (M is a 
normalization factor) 

P(M, ai ,...,a n ) =A A e -X 2 (M,a 1 ,...,a„) ) (31) 

where a\,...,a n is the relevant set of parameters and 

m(z;ai, ...,a n ) = M(M,H ) + 5log 10 (D L (z; ax, a n )) (3.2) 
is the apparent magnitude of the SNIa which depends on the luminosity distance 

Df(z; Oi, On) = (1 + z) f — - dZ ' . (3.3) 

J H(z';a 1 ,...,a n ) 

The quantity M is related to the absolute magnitude by 

M = M + 51og 10 (^)+25. (3.4) 

Moreover, for the cases of interest here, the Constitution dataset [9] with 397 data and the 
557 data points from the Union2 dataset [10], we have 

X (M,ai, ...,a n ) = 2^ 2 r~\ ' ( 

i=l a m obs\ Z i) 

where a m obs{zi) denotes the la errors of the SNIa data [15]. The parameters a\, ...,a n that 
minimize the \ 2 expression (3.5) are the most probable parameter values (the 'best fit') 
and the corresponding x 2 («i, a n ) = Xmin gi yes an indication of the quality of fit for the 
parametrization: the smaller xfnin 1S ' ^ ne better the parametrization. For our data analysis 
we follow the method developed in [15] which implies a marginalization over M and adapt 
it to the Constitution and Union2 data sets. The crucial point consists in an expansion for 
X 2 of equation (3.5) with respect to M of the form 

X 2 (ai,...,o n ) = A-2MB + M 2 C, (3.6) 



where 



397/557 , obs( , , h . 



A(ai,...,a n ) = 2^ 2 — r~\ ' ' 

m obs \ Z i) 



1 = 1 



. 39 ^ 57 K s (z i )- m tt (z i ;M = 0,a 1 ,..,fl n ) . . 

B( ai ,...,a n ) = ~ — i ( 3 - 8 ) 

i=l a m obs{ z i) 
397/557 

i=l a m °bs{Zi) 
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Equation (3.6) has a minimum for M = B/C at 

X 2 (ai, a n ) = A( ai , on) - g^igng . (3 . 10 ) 

Thus, instead of minimizing x 2 (M, a±, a n ) we can minimize x 2 (°i> ••> a n) which is inde- 
pendent of M. Obviously, xLn = Xmin- 

To be more specific, for the numerical analysis we use the data of table 1 (SALT) in [9] 
and of the Union2 data set [11]. To gauge our code, we reproduced in Fig. 2 the left panel 
of Figure 12 of [10] (SNIa + CMB) for a wCDM model with a constant EoS parameter oj 
without systematic errors and generalized it to the case of time-dependent EoS parameters. 



-0.2 




Figure 2. Reproduction (for SNIa+CMB) of the left panel of Fig. 12 in [10] for the wCDM model 
with two free parameters for SNIa+CMB. Our best-fit values (68.3% confidence level) are f2 m o = 
0.266±g:gi? and lo = -0.994±g;g^ compared with n m0 = 0.268lg$? and lu = -0.997^^55 in I 10 ]- 



3.1.2 BAO 

Acoustic oscillations at recombination give rise to a peak in the large-scale correlation function 
in luminous red galaxies [12, 16, 17]. This peak can be related to a distance scale 



Dv(zbao) 



ZBAO 



H(ZBAO) 

characterized by a dimensionless parameter 



ZBAO fa 



Mz b ao\ ai...a n ) = \/^m,oE(z B Ao) 



-1/3 







1 



H(z) 



ZBAO 



dz 



zbao Jo E(z; a\...a r , 



2/3 



(3.11) 



(3.12) 
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with zbao = 0.35 and E(z) = The observational value is A = 0.469 db 0.017 [12]. For 

Xbao we have 

2 / \ (A(z B Ao;ai...a n ) -0.469) 2 
X BAO (ai-a n ) = • (3-13) 

3.1.3 CMB 

The CMB shift parameter measures the displacement of the first acoustic peak of the CMB 
anisotropy spectrum with respect to the position this peak would have in a flat Einstein-de 
Sitter reference universe. For a flat universe it is given by [18, 19] 

r z is dz 

R(zi s \ax...a n ) = yO m / — r (3.14) 

J Q E[ K z;ai...a n ) 

with the last scattering redshift z\ s = 1090. The 7-years WMAP result for this parameter is 
R = 1.725 ± 0.018 [13]. For the Xcmb value we have 

= (3-15) 

For our analysis we combine the x 2 values from the three tests (3.10), (3.13) and (3.15) 
to minimize the total x\ tai 

Xtotal = XsNIa + XBAO + XCMB • (3.16) 

3.2 Analysis of the models 

For a non-interacting two-component model of the dark sector the consequences of a time- 
variable EoS parameter, based on the CPL parametrization and the Constitution data set, 
have been studied in [14]. For the Union2 set a similar analysis was performed in [21]. Our 
goal here is to extend this type of analysis to cosmological models with interactions in the 
dark sector. We shall rely on the CPL parametrization for the dark-energy EoS as well and 
use both the Constitution and the Union2 data sets. 

Our interacting models have 4 free parameters, namely the Hubble constant Hq, the 
fractional density of dark matter Q m o and the EoS parameters ojq and uj\. The baryon 
fraction Q^q is assumed to have the fixed value O^o = 0.042. As usual, Hq is parametrized 
by Ho = 100 h km sec _1 Mpc _1 . The parameter h will be marginalized according to (3.10). 
Hence, our models are left with the 3 free parameters £l m o, ujq and 

If baryons are taken into account, only the quantity £ldm0 = O m o — 0.042 interacts since 
baryons (with f^g = 0.042) are separately conserved. Generally, we find that the inclusion 
of baryons has only a very small influence on the curves. Therefore, baryons do not appear 
explicitly in our graphic representations. 

3.2.1 The ACDM model 

For later comparison we start with the best-fit data for the ACDM model in tables 1 and 2. 
These results are consistent with those of [25]. 
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Observation 


x 2 ■ 


^m0 


qo 


SNIa 


465.513 


0.289 


-0.566 


SNIa+BAO 


465.731 


0.282 


-0.577 


SNIa+CMB 


466.179 


0.278 


-0.583 


SNIa+BAO+CMB 


466.202 


0.276 


-0.585 



Table 1. Best-fit values based on the Constitution set for the ACDM model. 



Observation 


X 2 ■ 


^m0 




SNIa 


541.012 


0.269 


-0.596 


SNIa+BAO 


541.029 


0.271 


-0.593 


SNIa+CMB 


541.091 


0.266 


-0.601 


SNIa+BAO+CMB 


541.156 


0.268 


-0.598 



Table 2. Best-fit values based on the Union2 set for the ACDM model. 



3.2.2 The non-interacting case 

In a next step we consider the non-interacting case based on (2.10). The confidence level 
contours (la, 2a and 3a) for the Constitution set (short-dashed lines for SN only, long- 
dashed for a joint analysis SN+BAO and solid lines for the joint analysis SN+BAO+CMB) 
are shown in Fig. 3. The 2D contours are produced on the basis of 3 free parameters, i.e. 
A X 2 = X 2 ~ Xmin < 3.53(68%), 8.02(95%) and 14.02(99.73%). To obtain the contours in 
the ojq - ui plane (left panel), delta priors for O m o according to the best-fit values in Table 
3 were used. For the SN-only case, e.g., this means Q m o = 0.4519. The contours in the 
right panel rely on corresponding delta priors for uj\, i.e., u>\ = —11.227 for the SN-only case 
etc. These curves reproduce the results of [23]. The corresponding curves for the Union2 
sample are similar. The red arrow in the left panel characterizes the distance between the 
non-interacting model and the ACDM model (blue dot) for three free parameters if only the 
SNIa data are used. Under this condition, a value of Ax 2 ^ = 4.28 (cf. Table 5) corresponds 
to the 76.7% confidence region, equivalent to 1.31<r. Table 5 visualizes the relation between 
Ax 2 and the joint posterior probability P, given by [26] 

7 (i//2,Ax 2 /2) 

p=1 ~ nu/2) • < 3 - 17 > 

where v is the number of free parameters and ^ (v/2, Ax 2 /2) is the incomplete gamma 
function. 

The left panel of Fig. 4 shows the dependence of Xmin on ^mo for the range 0.20 < 
f^mo < 0.50 if only SNIa data (here, the Constitution set) are used. The minima correspond 
to much larger values of £l m o (0.452 for Constitution and 0.415 for Union2) than for the 
ACDM model (0.289 for Constitution and 0.269 for Union2). If combined, however, with 
BAO and CMB, these values reduce substantially, as shown in Table 3 (Constitution) and 
Table 4 (Union2). 

These tables also contain the best-fit values for ujq and uj\ for different data-set combina- 
tions and the corresponding values for the deceleration parameter. The right panel of Fig. 4 
displays the reconstruction of the deceleration parameter from the data. Note that using only 
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-2.0 -1.5 -1.0 -0.5 0.0 -1.5 -1.0 -0.5 0.0 



Figure 3. Non-interacting case: Contour plots (la, 2a and 3<r) in the ujq-uji plane (left panel) 
with delta priors for f2 m o from Table 3. Contour plots in the wo _ ^m0 plane (right panel) with 
corresponding delta priors for u\. The red arrow in the left panel characterizes the difference Ax 2 = 
X 2 (ACDM) — x 2 (Case No Int.) for the case of three free parameters (cf. Table 5). 




0.2 0.3 0.4 0.5 0.6 



"mo 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

z 

Figure 4. Left panel: Dependence of Xmin 011 ^mO for the Constitution set. Right panel: Dependence 
of the deceleration parameter on the redshift (with la contour lines) for two data-set combinations. 
For SNIa only, the best-fit model (O m0 = 0.4519, uj = -0.221, co 1 = -11.227 (cf. Table 3)) predicts 
a present decelerated expansion qo > 0. 



Observation 


y 2 • 


^m0 




UJl 


Qo 


SNIa 


461.231 


0.452 


-0.221 


-11.227 


0.318 


SNIa+BAO 


465.425 


0.281 


-0.905 


-0.497 


-0.475 


SNIa+CMB 


465.557 


0.270 


-0.935 


-0.119 


-0.524 


SNIa+BAO+CMB 


465.606 


0.274 


-0.965 


0.015 


-0.550 



Table 3. Best-fit values based on the Constitution set for the non-interacting model with the Hubble 
rate (2.10). 
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Observation 


v 2 . 


^m0 


w 




qo 


SNIa 


539.878 


0.415 


-0.886 


-5.108 


-0.278 


SNIa+BAO 


540.988 


0.274 


-1.007 


-0.039 


-0.598 


SNIa+CMB 


541.028 


0.266 


-1.010 


0.103 


-0.612 


SNIa+BAO+CMB 


541.070 


0.269 


-1.030 


0.187 


-0.630 



Table 4. Best-fit values based on the Union2 set for the non-interacting model with the Hubble rate 
(2.10). 



the SNIa data from the Constitution set, the best-fit model (f2 m o = 0.4519, uq = —0.221, 
uj\ = —11.227 (cf. Table 3)) predicts a decelerated expansion q$ > (cf. [22]). This is differ- 
ent from the result qo < based on the Union2 sample (cf. Table 4) All the results for the 
non-interacting model, based on the Constitution set, are consistent with those of [14] and 
[23], where the latter reference also discusses the tensions between the different SNIa data 
sets. 



Observation 


No Int. 


ACDM 




P 


SNIa 


461.231 


465.513 


4.282 


76.7% (1.31a) 


SNIa+BAO 


465.425 


465.731 


0.306 


4.1% (~0.06a) 


SNIa+CMB 


465.557 


466.179 


0.622 


10.9% (~0.16o-) 


SNIa+BAO+CMB 


465.606 


466.202 


0.596 


10.3% (~0.15o-) 



Table 5. Minimum \ 2 values for the non- interacting and the ACDM models and the differences 
— X 2 (ACDM) — x 2 (Case No Int.). The right column provides the corresponding probability 
values for the case of three free parameters. For comparison we recall that A\ 2 = 3.53 corresponds 
to P = 68.3%(1ct), Ax 2 = 8.02 corresponds to P = 95.4%(2ct) and A X 2 = 14.2 corresponds to 
P = 99.73%(3ct) [26]. 

3.2.3 The case £ = 1 

The results for this model are shown in Figs. 5 and 6, the corresponding best-fit values are 
given in tables 6 and 7. Different from the non-interacting case, the minimum value of x 2 f° r 
the SNIa samples is practically independent of £l m o- But there is a slight tendency to f2 m o ~ 
(left panel of Fig. 6). This degeneracy is also seen in the solid lines of the right panel of Fig . 5 
where only SN data (Constitution) have been used. It shows elongated (with respect to fi m o) 
contour plots and has more negative values for uq than the non-interacting case. For the 
SNIa analysis we included in tables 6 and 7 also the results for a prior J7 m o = 0.289. There 
is a tendency to larger values of u\ compared with the non-interacting case. The left panel 
of Fig. 5 reveals that both data set combinations indicate oj\ ~ 1 while ujq ~ — 1 for the joint 
analysis, something different from the non-interacting case as well. The insert magnifies the 
relevant region. The deceleration parameter behaves almost in the same way for the different 
data combination as can be seen in the right panel of Fig. 6. 

3.2.4 The case £ = 3 

This model is characterized by Figs. 7 and 8 as well as by tables 8 and 9. Again, the minimum 
value of x 2 f° r t ne SNIa samples varies only slowly with Q m o with a slight tendency toward 



- 11 - 



-1.0 
io 



-0.5 0.0 



-1.5 -1.0 -0.5 0.0 

li<0 



Figure 5. Case £ = 1: Contour plots (la, 2a and 3a) in the ujo-uji plane (left panel) with delta priors 
for fi m o from Table 6. Contour plots in the ojQ-Q m Q plane (right panel) with delta priors for wi from 
Table 6. 
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Figure 6. Case £ = 1. Left panel: dependence of Xmin on ^mo- Notice that Xmin remains practically 
constant over the entire range of fl m o- Right panel: dependence of the deceleration parameter on the 
rcdshift (with la contour lines) for two data-set combinations. 
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SNIa 


465.906 


0.000 


-0.753 


1.051 


-0.629 


(Prior Q m0 = 0.289) 


466.022 


0.289 


-1.054 


1.167 


-0.624 


SNIa+BAO 


466.018 


0.276 


-1.035 


1.159 


-0.624 


SNIa+CMB 


466.015 


0.266 


-1.018 


1.135 


-0.621 


SNIa+BAO+CMB 


466.020 


0.275 


-1.039 


1.191 


-0.630 



Table 6. Best-ht values based on the Constitution set for the case £, = 1 with the Hubble rate (2.15). 



O m o = 1 (left panel of Fig. 8). Therefore, the first rows in tables 8 and 9 remain empty since 
no reliable entries are available. For the same reason there do not appear confidence contours 
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On 


SNIa 




541.223 


0.000 


-0.779 


1.045 


-0.669 


(Prior Q m0 = 


= 0.269) 


541.275 


0.269 


-1.063 


1.171 


-0.665 


SNIa+BAO 




541.275 


0.271 


-1.066 


1.173 


-0.665 


SNIa+CMB 




541.272 


0.252 


-1.037 


1.153 


-0.664 


SNIa+BAO+ 


-CMB 


541.300 


0.270 


-1.081 


1.269 


-0.683 



Table 7. Best-fit values based on the Union2 set for the case £ = 1 with the Hubble rate (2.15). 



for the SNIa samples only in Fig. 7. In the second rows of tables 8 and 9 we include the 
parameter values for the priors f2 m o = 0.289 and Q, m o = 0.269, respectively. The preferred 
wi-values are closer to the non-interacting case than to the case £ = 1 (left panel of Fig. 7) . 
Again, an insert magnifies the relevant region. Also the elongation in the right panel of Fig. 7 
is reduced compared to its ^ = 1 counterpart and, as in the non-interacting case, ujq > — 1 is 
preferred. 

At this place a comment on the redshift dependence of the deceleration parameter in 
the right panels of Figs. 4, 6 and 8 is in order. According to the joint analysis, q{z) decreases 
with decreasing z in all the cases. This contrasts with the results found in [22], according 
to which there exists a minimum in q{z) which would imply a slow down of the accelerated 
expansion and, possibly, a transition back to decelerated expansion. However, as Fig. 4 also 
shows, considering only the SN data, the non-interacting model shows the same slowing- 
down behavior as that found in [22]. From the Union2 data alone, a tendency that the 
cosmic acceleration will slow down was also reported in [20] . For a combination of SNIa with 
the BAO data D v {z = 0.35)/ D v (z = 0.20) this result was confirmed in [22]. For an analysis 
for D v (z = 0.35), however, we were unable to reproduce this minimum. Also, if CMB data 
are additionally included, the mentioned tendency for q(z) to grow with decreasing redshift 
close to z = disappears. A detailed analysis of this type of tensions has recently been 
performed in [21]. A discussion on different kinds of tensions, including those between SNIa 
and CMB data, can be found, e.g., in [23, 24]. 
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** 


** 


** 


** 


** 


(Prior Q m0 = 0.289) 


465.482 


0.289 


-0.971 


-0.187 


-0.345 


SNIa+BAO 


465.495 


0.277 


-0.958 


-0.126 


-0.414 


SNIa+CMB 


465.497 


0.277 


-0.961 


-0.093 


-0.451 


SNIa+BAO+CMB 


465.497 


0.277 


-0.963 


-0.088 


-0.457 



Table 8. Best-fit values based on the Constitution set for the case £ = 3 with the Hubble rate (2.19). 



3.2.5 Comparing the models 

The different models may be properly compared among themselves and also with the ACDM 
reference model, using appropriate statistical criteria. Two options are the already mentioned 
AIC and BIC criteria, which allow us to compare models with a different number of degrees 
of freedom. The AIC criterion uses the formula AIC = Xmm + ^/c [27], where k is the number 
of degrees of freedom; the BIC criterion [28] is based on the expression BIC = Xmin+^k In N, 
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Figure 8. Case £ = 3. Left panel: dependence of Xmin 011 ^mO- Notice that also here Xmin remains 
practically independent of f2 m o- Right panel: dependence of the deceleration parameter on the redshift 
(with lcr contour lines) for the combination SNIa + BAO + CMB. 



Observation 


X 1 ■ 


^m0 








SNIa 


** 


** 


** 


** 


** 


(Prior n m0 = 0.269) 


540.996 


0.269 


-1.017 


0.105 
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SNIa+BAO 


540.994 


0.272 


-1.020 


0.083 


-0.694 


SNIa+CMB 


540.994 


0.272 


-1.022 


0.101 


-0.715 


SNIa+BAO+CMB 


540.997 


0.272 


-1.018 


0.092 


-0.702 



Table 9. Best-fit values based on the Union2 set for the case £ = 3 with the Hubble rate (2.19). 



where iV is the number of observational points. The smaller the resulting numbers in both 
expressions, the higher the quality of the corresponding model. It is convenient to classify 
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a model with respect to the differences AAIC and ABIC between its AIC and BIC values, 
respectively, and the corresponding values for a reference model. This establishes a scale 
which allows for a ranking of different models according to the magnitude of their differences 
AAIC and ABIC [29, 30]. The smaller the difference to the lowest AIC or BIC values, 
here those of the ACDM model, the better the model. For differences less than 2, there 
is strong support for the model under consideration. If AAIC(ABIC) < 6 the model is 
still weakly supported. Models with AAIC (ABIC) > 10 should be considered as strongly 
disfavored. In tables 10 and 11 we summarize the results for the investigated models and 
assess them according to the AIC and BIC criteria. Notice that all the Xmin values for the 
competing models are smaller than that of the ACDM model. But the mentioned criteria 
penalize the introduction of additional parameters and reverse the ranking. By inspection, 
it follows, that, using the AIC criterion, the non- interacting model and the £ = 3 model 
are still weakly supported. Applying, however, the BIC criterion, all these models are ruled 
out. This kind of contradiction in using different evaluation criteria is well known in the 
literature, see, e.g., [31]. Graphical summaries of our analysis is given in Figs. 9 and 10. 
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q = -0.630 


q Q = -0.457 


X 2 ■ 


466.202 
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3.404 


3.818 
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Table 10. Summary of the analysis for the Constitution data set. 
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n m0 = 0.272 
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q = -0.630 


qo = -0.683 


q = -0.702 


Arm?i 


541.156 
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3 


3 
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12.559 


12.789 


12.486 


AAIC 





3.914 


4.144 


3.841 



Table 11. Summary of the analysis for the Union2 data set. 



Finally, we consider the direction of the energy transfer between the dark components. 
According to (2.7), the sign of Q for the current Universe depends on the sign of £/3 + coq. 
As to be seen in Fig. 11, our interacting model with £ = 1 corresponds to Q > for both 
combinations of data sets. The case £ = 3 lies, as expected, very close to the non-interacting 
model. The Constitution data set yields Q > while the Union2 data set slightly favors an 
energy transfer from dark matter to dark energy, i.e, Q < 0. 
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4 Summary 



Cosmological models in which an interactions between dark matter and dark energy is ad- 
mitted, give rise to a richer cosmological dynamics than non-interacting models, albeit at 
the expense of an additional parameter. The introduction of interacting models is largely 
motivated by the possibility to address the coincidence problem. In this paper we have in- 
vestigated a class of interactions that result in a power-law behavior p m /p x oc of the 
ratio of the energy densities of dark matter and dark energy. Generalizing previous work, 
we admitted a time-varying EoS parameter of the dark-energy component within the CPL 
parametrization. We found analytic solutions for the cases £ = 1 and £ = 3. The former is of 
interest with respect to an alleviation of the coincidence problem, the latter primarily to test 
a potential time variation of the EoS parameter. With the help of a Bayesian statistical anal- 
ysis we tested the resulting dynamics against the SNIa data of the Constitution and Union2 
samples. We included also information from BAO and CMB shift data and compared the 
interacting models among themselves and with a non-interacting model. The Xmin v& l ues 
for all the competing models turned out to be smaller than the Xmin varue f° r the ACDM 
model. But according to both the AIC and the BIC criteria which penalize the introduction 
of additional parameters, the ACDM model remains the preferred choice. 

Our study was restricted to the homogeneous and isotropic background dynamics. A 
more complete analysis requires to investigate the implications for structure formation as 
well. We believe, that the analytic solutions for the Hubble rates, found in this paper, will be 
helpful to calculate the matter power spectrum and the impact on the integrated Sachs- Wolfe 
effect in future work. 
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Figure 9. Graphical summary I of the analysis. 
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Figure 10. Graphical summary II of the analysis. 
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Figure 11. Direction of the energy transfer at the present epoch. The straight line divides the regions 
Q > and Q < 0. The points denoted by (1) correspond to the set Constitution+BAO+CMB, the 
points denoted by (2) correspond to Union2+BAO+CMB. The inserts magnify the la regions. 



